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Clues to the shape of our Universe can be found by searching the CMB for matching circles of 
temperature patterns. A full sky search of the CMB, mapped extremely accurately by NASA's 
WMAP satellite, returned no detection of such matching circles and placed a lower bound on the 
size of the Universe at 24 Gpc. This lower bound can be extended by optimally filtering the WMAP 
power spectrum. More stringent bounds can be placed on specific candidate topologies by using 
a a combination statistic. We use optimal filtering and the combination statistic to rule out the 
infamous "soccer ball universe" model. 



I. INTRODUCTION 

What is the shape of space? While this question 
may have once seemed more philosophical than scientific, 
modern cosmology has the chance to answer it using the 
oldest observable light in the Universe, the Cosmic Mi- 
crowave Background radiation (CMB). NASA's Wilkin- 
son Microwave Anisotropy Probe (WMAP) has made a 
detailed map of the CMB sky which has been used to pro- 
vide answers to many age-old questions about the nature 
of the Universe |2| ■ 

While it is certainly possible that the Universe extends 
infinitely in each spatial direction, many physicists and 
philosophers are uncomfortable with the notion of a uni- 
verse that is infinite in extent. It is possible instead that 
our three dimensional Universe has a finite volume with- 
out having an edge, just as the two dimensional surface 
of the Earth is finite but has no edge. In such a universe, 
it is possible that a straight path in one direction could 
eventually lead back to where it started. For a short 
enough closed path, we expect to be able to detect an 
observational signature revealing the specific topology of 
our Universe [l|. 



II. GEOMETRY AND TOPOLOGY 

An important question answered by the WMAP mis- 
sion is that of the curvature of space. The matter and 
energy density of the Universe indicate that space is very 
nearly flat. The WMAP data point to a universe with a 
total energy density within 2% of critical j|i|. This means 
that even if space in not quite flat, the radius of curva- 
ture of the Universe is at least of order the size of the 
observable Universe, and space can be considered to be 
nearly flat. 

The WMAP sky also provides clues about the topol- 



ogy of the Universe. Observable topologies would have a 
characteristic fingerprint in the CMB sky. The set of pos- 
sible topologies for our Universe is determined by the cur- 
vature of space. In a flat Universe, the allowable topolo- 
gies are restricted to a set of eighteen possibilities. It has 
been shown that a nearly flat Universe would have an ob- 
servational fingerprint very similar to that of an exactly 
flat Universe [3,l3- Cornish, Spergel, and Starkman [l| 
have described the CMB signature revealing the shape of 
space as "circles in the sky" . 



III. CIRCLES IN THE SKY 

Circles in the CMB revealing the topology of the Uni- 
verse could only be observed for a universe that admits 
sufficiently short closed paths. If our Universe is mul- 
tiply connected and the fundamental domain fits inside 
the observable Universe, we expect to see the character- 
istic signature of the topology in the CMB. In a Universe 
with a non-trivial topology, there exist multiple copies of 
objects, one in each copy of the fundamental domain. A 
copy of our own galaxy, solar system, and planet could 
possibly be observed many light-years away. Surrounding 
each copy of Earth would be a copy of the last scatter- 
ing surface - a 2-sphere at redshift z ~ 1100. Since in a 
homogeneous and isotropic geometry the intersection of 
two 2-spheres is a circle, and these circles in the CMB are 
physically the same place in space, to the extent that the 
CMB measures the local gravitational potential, density 
and temperature on the last scattering surface, we expect 
to find matching circles of temperature patterns in the 
CMB when looking in two different directions in the sky 
(Figure P). 

The search for these matching circles of temperature 
patterns in the WMAP data is a huge task. The resolu- 
tion of the WMAP satellite provides over three million 
independent data points on the sky. Each of these points 




FIG. 1: In a multiply connected universe, all of space must 
be tiled with copies of the fundamental domain. Here the 
fundamental domain is a cube, and associated with each do- 
main is a copy of the last scattering surface. The intersections 
are physically the same place in space, producing matching 
circles of temperature patterns in two different directions as 
observed from the Earth. 



Here Tim denotes the m*^ harmonic of the tempera- 
ture pattern around the i^^ circle. The parameters a 
and (3 denote the angular radius and the phase offset 
between the circles, respectively. The S statistic is nor- 
malized with a perfect match yielding an S value of 1. 
The comparison is done in Fourier space, in order to in- 
clude appropriate weighting of different angular scales. 
The m weighting factor on circles plays the same role 
as the usual i{2i + 1) weighting that appears in the full 
angular power spectrum: these factors take into account 
the numbers of degrees of freedom per mode [3 • If the m 
weighting is neglected large angular scales dominate and 
completely drown out small-scale fluctuations. A limited 
search for matching circles in the CMB sky has also been 
performed in Ref. ^. There, the matching statistic of [l| 
was used; one that does not include the appropriate m 
weighting factor: 



can be considered a circle center, with the circle radius 
ranging between and 90°. Once the temperature pat- 
tern around a given circle is found, it must be compared 
to every other circle in the sky with the same angular ra- 
dius, and the matching statistic calculated for the pair. 
The phasing of each circle pair must also be considered, 
ranging from to 360°, since the start of the compari- 
son can begin anywhere along the second circle. This full 
search is time consuming, and is still in progress. 

A more efficient search can be employed when consid- 
ering the WMAP result of nearly zero spatial curvature 
for the Universe. In a flat universe, all eighteen possi- 
ble topologies produce some circle pairs that are back- 
to-back (though the smallest circles may not be back- 
to-back). That is, their circle centers are 180° apart. 
It can also be shown that for a large radius of curva- 
ture compared to the radius of the last scattering surface, 
non-trivial topologies with a positive curvature produce 
nearly back-to-back circle pairs |^ |^ . A shorter search 
can thus be used on the WMAP data, where each circle is 
compared only to its back-to-back pair. Similarly, an al- 
most back-to-back pair search can be implemented, look- 
ing only for matching circles between 170° and 180° from 
the first circle center. This search has been completed, 
and no statistically significant matches were found Q. 



IV. THE S STATISTIC 

A matching statistic is needed to define the correla- 
tion between the temperature patterns around two circles 
with the same angular radius. In Ref. '4'|, the S statistic 
for circle comparison is defined, 
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Here the brackets denote integration along the circle. 
This alternative definition of the 5* statistic can be calcu- 
lated in position space, but the consequence of disregard- 
ing the weighting factor is that large scale temperature 
variations swamp any small scale features. The impor- 
tance of the m weighting is apparent in Figure |21 
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FIG. 2; The Fourier coefficients of 6T, as a function of mul- 
tipole moment m, around a circle with no m weighting (top) 
and including m weighting (bottom). 



V. FALSE POSITIVE CALCULATION 

An essential aspect of our search is the calculation of 
the threshold for a significant match. We expect this 



false positive line to be higher for circles with smaller an- 
gular radius, since these circles have fewer pixels used for 
comparison and thus have a greater chance for random 
matching. It has been shown in Ref. ^ that the level for 
a positive detection can be calculated using 
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Here iVgcarch is the size of the search space, counting the 
number of circle comparisons made in the search at each 
angular radius, a. The false positive threshold is sensitive 
to A'soarchj increasing for larger searches. The expected 
value of S for random skies can be estimated from first 
principles, but it is better to numerically measure the ex- 
pectation value from the data. The expectation value of 

S, (5*^) ^ , and the size of the search space, A^search, follow 
from the implementation of the specific search being per- 
formed. When looking at the results from our searches, 
we plot a false positive threshold for which we expect 
fewer than 1 in 100 random skies to produce an S value 
above the line. 

Calculating the correct value of iVscarch involves deter- 
mining the number of completely independent possible 
searches. For example, a sky with 2° smoothing has less 
independent pixels than a sky with 1° smoothing and 
thus has fewer possible independent circle pairs to be 
searched. The full iVgcarch value for a search should take 
into account the number of independent circle radii (a) 
searched at a given sky resolution. This was not included 
in Ref. [j|, but we include it here by multiplying the 
-^search valuc at each a calculated by the search code by a 
factor equal to the number of independent radii searched. 
The false positive level must be considered for any claim 
of detection, and is used to place the lower bound on the 
size of the Universe for no circle pair detections. 



VI. THE WMAP BOUND 

The full sky search for back-to-back and almost back- 
to-back matching circles was performed in Ref. Q, and 
it is apparent in Figures 3 and 4 of that paper that no 
statistically significant matches were found. 

A simulation of the CMB has been produced |j] that 
includes all relevant physics, using model parameters that 
give a good match to the real CMB power spectrum, but 
with a non-trivial topology built in. This 3-torus universe 
with a fundamental domain that fits inside the sphere of 
the CMB is used to test the search codes. The results of 
the back-to-back search on the simulated 3-torus universe 
indeed show peaks in the matching statistic, indicating 
matched circle pairs (Ref. |j|. Figure 1). The fact that 
the results of the same back-to-back and nearly back-to- 
back searches of the WMAP data show no such peaks can 
be used to place a lower limit on the size of our Universe. 
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FIG. 3: Contributions to the power spectrum of the CMB. 
The effective temperature around circles can be correlated due 
to topology, but is degraded by ISW effects. The acoustic ve- 
locity also degrades correlation for small circles. (Reproduced 
with permission from Nature, Hu, Sugiyama, & Silk (1997)) 



The intersection of the false positive threshold and the 
matching statistic peaks found for a simulated 3-torus 
universe defines this lower limit, since it indicates that 
the smallest circle pairs we could expect to detect would 
be approximately 20°. 

For decreasing a, a downward trend can be seen in the 
3-torus simulation results for the height of the peak indi- 
cating a match. Photons from the last scattering surface 
travel toward Earth, carrying with them an imprint of 
the conditions at their point of origin. It is this effec- 
tive temperature of the CMB photons that should match 
along a circle. However, two matching photons travel- 
ing different paths toward Earth will encounter different 
line of sight effects, degrading the temperature match. 
The combined line of sight effects on CMB photons can 
be understood by considering the physical contributions 
to the total power spectrum of the CMB (see Figure |3J|. 
The Integrated Sachs- Wolfe effect (ISW) degrades any 
possible match, while the contribution by the acoustic 
velocity term depends on the size of the circle pair. The 
velocity term is given by n • w , where n is a unit vec- 
tor along the line of sight and v is the velocity vector of 
the plasma at the surface of last scatter. For large cir- 
cles, this term is correlated, while for small circles this 
term becomes increasingly anti-correlated as the circles 
get smaller, (see Figure QJ. This velocity term thus de- 
grades the matching statistic value for small circles and 
explains the downward trend of peaks for small a. 

The null result of the search for matching circles in 
the CMB indicates that the fundamental domain of the 
Universe must be at least on order the size of the surface 
of last scatter. The fact that the false positive line in- 
tersects the maximum peaks expected for the matching 
statistic at a = 20° means that the fundamental domain 




FIG. 4: Here a circle represents a cross section of the 2-sphere 
of the CMB. The point of intersection thus represents a point 
along a circle in the sky, viewed from two different copies of 
the Earth. It can be seen that small circles have an almost 
anti-parallel velocity contribution, h ■ v (top), while large cir- 
cles have an almost parallel velocity term (bottom). 



must be big enough such that only circles smaller than 
20° could be produced by intersections of copies of the 
CMB. This, along with the best fit values for the other 
cosmological parameters, places a lower bound on the size 
of the Universe at 24 Gpc. 



VII. EXTENDING THE BOUND 

There are a couple of techniques that can be applied 
in an attempt to extend this bound beyond 24 Gpc (or 
possibly detect the topology of the Universe, lurking just 
below the current false positive threshold). Since we are 
interested in raising S values for small circles where it 
dips below the detection line, we can filter out from the 
CMB power spectrum the terms that degrade matches for 
small circles. This requires filtering out both the ISW 
effect on large angular scales and the acoustic velocity 
term on smaller angular scales. We thus multiply the 
WMAP power spectrum by a filter that includes only 
the regions of the power spectrum where the effective 
temperature dominates all other terms. 

The basic shape of our filter comes from the physical 
theory described in Figure |31 We then search for the 
best filter on the 3-torus simulation, varying the location, 
width, and steepness of the filter windows. Our best 
simulated filtered power spectrum is seen in Figure |S1 

The filtering indeed raises the S values for small circles 
in our simulated 3-torus universe and lowers the S values 
for large circles, as expected when filtering out the veloc- 
ity term that contributes to raising S for large angular 
radius pairs (Figure EJ. 

In general, the false positive line depends on the filter 
being used. A filter that keeps only the low I power, 
for instance, would result in a higher false positive line 
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FIG. 5: The unfiltered simulated power spectrum (dotted 
line) and the power spectrum with regions dominated by ei- 
ther ISW or acoustic velocity effects filtered out (solid line). 
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FIG. 6: The dashed lines give the S values for the unfiltered 
simulation, while the solid lines give the S values when the 
power spectrum filter is included. Top: The average values 
of the S statistic as a function of angular radius. Bottom: 
The maximum values of S as a function of angular radius. 
Smax increases by ~ 10% for small circles with filtering, but 
decreases for large circles, as expected. Plotted with the cal- 
culated false positive line. 



since there are essentially less significant pixels without 
the high I values, and thus a greater chance of a false 
match. Our power spectrum filter keeps most of the high 
I values, especially near the first acoustic peak, and thus 
does not significantly affect the false positive line. 

Another choice of filter is the use of m weighting, as 
discussed in Section Hvl The S statistic calculation ne- 
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FIG. 7; A back-to-back circle search is performed on our 
3-torus simulation with 1° degree smoothing. The dashed 
line shows the results for a search with the matching statistic 
calculated without the necessary m weighting, while the solid 
line indicates a search including m weighting. The threshold 
line for detection for each search is also plotted. 



glecting m weighting results in higher S values for circle 
pair comparisons at all angular radii, raises the false pos- 
itive level, and degrades significant peaks. These effects 
are illustrated in Figure [7| showing the results of a back- 
to-back circle search of a 3-torus CMB simulation when 
m weighting is neglected (dashed line) and included (solid 
line). It can be seen that the search without m weight- 
ing emphasizes large angular scales, and thus only finds 
matching circles with a > 45°. 

Several recent papers [llLll3 . ll^ have called into ques- 
tion the efficacy of the matched circle test. In each in- 
stance the negative conclusions can be traced to the use 
of the unweighted S statistic. As our searches of real- 
istic simulated skies have shown, the properly weighted 
matching statistic provides a powerful tool for probing 
the topology of the universe. 



VIII. COMBINATION STATISTIC 

Another technique for probing below the limit of 20° 
circle pairs is employing a directed search for a specific 
topology. Such a search can increase the S value for 
small circles by stringing together several sets of circles in 
an expected configuration to create an effectively larger 
radius circle with the increased number of pixels. For 
a cubic 3-torus universe, for example, we expect to find 
at least three sets of circle pairs corresponding to the 
matching faces of the fundamental domain and we can 
predict the relative orientations of each matching set. 

The combination search for this example would con- 
sist of choosing a circle center and circle radius and com- 
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FIG. 8: The combination statistic calculated for our 3-torus 
simulation using the optimal filter (solid line). The upper 
dashed line gives the threshold for less than 1 in 100 random 
skies producing an S value above the line. 



paring the pixels around this circle to the circle exactly 
opposite it in the sky. The search then chooses an axis 
for a second set of back-to-back circles with the same an- 
gular radius. The location of the final circle pair is fixed 
once the axis has been chosen. The pixels along all three 
pairs are strung together to calculate the total matching 
statistic for the whole set. Since there are more pixels 
for the whole set than for one pair of circles, there is a 
smaller chance for a random match. 

This combination statistic can be tested using our sim- 
ulated 3-torus universe. There are matching circles with 
angular radius of about 16° in our simulation. These 
pairs cannot be detected in the single pair back-to-back 
circle search because the false positive threshold is above 
the peak of their matching statistic. 

Our 3-torus simulation has a fundamental domain such 
that several copies fit within the surface of last scatter. 
There are thus more than three sets of circles at 16°, 
and our combination statistic is calculated by stringing 
together the pixels along four sets of circles. 

It can be seen in Figure El that the search for the com- 
bination statistic in our simulated 3-torus universe raises 
the S value above the false positive line, and we can 
thus make a statistically significant detection. The 3- 
torus combination statistic can be applied to the real 
WMAP data to rule out (or detect) sets of matching cir- 
cles smaller than the 20° cutoff of the single circle pair 
search. 

It is crucial to calculate the false positive line cor- 
rectly for the combination search. While the combination 
statistic lowers the false positive level for small circles by 
creating effectively larger circles, it must be remembered 
that the combination statistic requires a larger search 
than looking for back-to-back circles. There is an extra 
degree of freedom for the orientation of the fundamen- 
tal domain once the first set of circles is identified. The 
false positive level for the combination search can be es- 
timated from the standard pair-wise search as follows. 




FIG. 9: The Poincare dodecahedron is a 12-sided polygon 
with identical faces. This shape can tile space with no gaps 
only in a positively curved universe. The face identifications 
include a twist of ±36° . 



The RMS value for fi-eombo is given by (S^^^bo)'^' = 
(Sp^ij.)^/^/?!^/^, where n is the number of circle pairs 
used in the combination statistic. The number of circle 
comparisons in the combination search, A'j!g°™j^j[' scales as 
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, where d is the number of degrees of freedom 
in the pair-wise search. Combining these observations in 
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We can also estimate the combination statistic level for 
matching circles in terms of the pairwise values. The 
combination statistic for matching circles, 5'™^^^, is the 
weighted average of the pair-wise statistic for matching 
circles, 5'™^''^'^, and since the total power around each 

circle is roughly equal, we have S™f^^ ~ .^^fj'^'^ where 
the bar denotes the average. These estimates for the false 
positive and match level for the combination searches 
were found to be accurate predictors of the numerical 
results. Extrapolating from the values of 5™^^'^'' seen in 
the upper panel of FigureEland the values of 5'™f*'^'^ seen 
in the lower panel of Figure[7|suggests that a combination 
search using 4 of more circle pairs should be able to detect 
circles with radii as small as 5°. 



IX. POINCARE DODECAHEDRAL SPACE 

In a universe with positive curvature, the faces of a 
dodecahedron can be identified to form what is known 
as Poincare Dodecahedral Space (the "soccer ball uni- 
verse"). The face identifications for a dodecahedron in- 
volve a twist of ±36° to insure the matching up of the ver- 
tices of the face pair (Figure O. It has been claimed that 
the power spectrum expected for a universe with such a 
topology closely matches the power spectrum found for 
our Universe by WMAP ^] . 



X. CLAIM OF A NON-TRIVIAL TOPOLOGY 

A claim has been made that six sets of matching circles 
in the CMB have been found, indicating that the Uni- 
verse has the topology of the Poincare Dodecahedron J3] ■ 
The claim, however, is not accompanied by a statistical 
analysis of the results, and cannot be assessed until this 
is completed. 




FIG. 10: Our combination statistic plotted for circle compar- 
isons with a —36° twist (solid line) and a ±36° twist (dashed 
line) . Including 2° smoothing, cutting out the galaxy contam- 
ination, and neglecting the m weighting factor reproduces the 
peak found in Ref. j^. 

In Ref. ^, the combination statistic for the expected 
six sets of matching circles of a Poincare dodecahedron 
has been computed for a range of circle radii a with a 
peak appearing for circles around 11°. The peak only ap- 
pears when the —36° twist is used when identifying faces 
of the dodecahedron (c.f. Ref. S|, Figure 4). This peak is 
then compared to the same search using the ±36° twist 
for face identifications (c.f. Ref. Q], Figure 5), where 
there appears to be no peak. It is expected that the 
peak would appear for only one twist value if this is an 
indication of the physical topology of the Universe. 

Figure ^| shows that we were able to find the same 
peak in the matching statistic around 11° using a com- 
bination statistic. In order to reproduce the peak, we 
used the S statistic given in Equation ^, a galaxy cut, 
and the 2° smoothing used in Ref. 8]. The locations of 
the six pairs in the CMB sky, shown in Figure [TTl also 
matches those found in Ref. l^j. 

For proper analysis of the results, the false positive 
line must be calculated correctly. The results from the 
±36° twist search do not give the actual false positive 
line for detection. To get a better idea of the expected 
peak levels for a universe with a trivial topology, the 
WMAP data can be randomly scrambled. This preserves 
the power spectrum and distribution of temperature val- 
ues while destroying any matching circles that would give 
an indication of topology. The same search can then be 
performed on many of these scrambled skies and the re- 
sults plotted together to show the expected random peak 




FIG. 11: The locations of the six sets of circle pairs in the 
CMB sky identified in Ref. 8ij and found by our combination 
statistic search. 




FIG. 12: Top: The combination search for the dodecahedral 
topology performed with several different twist values. The 
expected level for S is plotted in bold. Bottom: The possible 
detection with the —36° twist (solid line) is plotted with the 
appropriate expected level for S (lower dashed line) and the 
threshold for fewer than 1 in 100 random skies producing an 
S value above the line (upper dotted line) 



levels for a trivial topofogy. In the case of the position 
space search used in Ref. Q, the WMAP sky cannot 
be randomly scrambled while preserving both the power 
spectrum and the galaxy cut. To get a reliable measure 
for the expected S values for the combination search we 
instead use various twists other than ±36° when identi- 
fying circle pairs. The results from the search on these 
random twists is plotted along with the possible detection 
in Figure El 

It can be seen in Ref. Q, Figures 13 through 18, that 
two sets of circles include points that run through the 
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FIG. 13: The distribution of S values calculated at a = 11° 
(top) and a = 20° (bottom). Both have a non-zero mean due 
to a DC component of the temperature values around circle 
pairs. 



plane of our galaxy. The CMB data for this portion of 
the sky is contaminated by the radiation from our galaxy. 
While it is possible to filter out much of this contamina- 
tion, a more conservative approach is to perform a galaxy 
cut. The galaxy cut is easily done in position space, but 
is not quite as easily included in a Fourier space search. 
Including the galaxy cut raises the peak for matching 
circles, since it cuts out uncorrelated points, but it also 
raises the false positive level, since there are less data 
points to compare. It can be seen in Figure [HI that the 
unweighted Fourier and position space searches produce 
similar peaks and that the galaxy cut in position space 
raises the peak as expected. While Parceval's theorem 
states that the position space search with no galaxy cut 
and the Fourier search with no m weighting should pro- 
duce exactly the same result, numerical error leads to a 
slight difference in the results. 

We expect that the S values calculated for all of the 
sets of circle pairs compared at a give a should form a 
Gaussian distribution. Even if there is an outlier value 
of Smax indicating the topology of the Universe, this one 
value should not significantly skew the distribution of all 
of the non-detections at the same angular radius. Two 
distributions of calculated S values for the —36° twist 
search can be seen in Figure El The non-zero mean of 
the distributions indicates an average value for the dis- 
tribution above zero, and thus a DC contribution to the 
temperature values around each circle pair. This con- 
tribution can be subtracted to show the true compari- 
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FIG. 14: The search in Fourier space with no m weight- 
ing (solid line) matches the search in position space with no 
galaxy cut and no m weighting (dashed line). Including the 
galaxy cut raises the peak (dotted line). 



FIG. 16; Our results for the Poincare Dodecahedron combi- 
nation search with optimal filtering. The dashed line indicates 
the detection threshold. 
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FIG. 15: Removing the DC contribution (average value) 
when comparing circle pairs shifts the distribution of S values 
toward zero mean and indicates no peak at 11°. 



son of temperature variation. Simply removing the DC 
contribution reduces the peak in the S values of the do- 
decahedron combination search (Figure |15I top panel). 
The distribution of S values is shifted toward zero mean 
by explicitly removing the average temperature values, 
but the mean is not exactly zero. This is due in part 
to numerical error. The S statistic calculated using m 
weighting (Equation (Q) automatically takes care of this 
average offset. 

Figure ITCl shows our final analysis, using the Poincare 



combination search with the appropriate m weighting 
and plotted with the threshold line for detection. The 
appropriate analysis yields no significant peak at 11°. It 
is thus clear that when the the proper matching statistic 
is used, and the expected level of false positives is calcu- 
lated, that the detection claimed in Ref.g^J does not hold 
up. It seems that the observed peak is most likely due 
to overemphasis of large angular scales by neglecting m 
weighting, resulting in searching for matching circles in a 
sky with large patches of constant temperature. Such a 
sky would indeed produce randomly matched small cir- 
cles. 

It should also be noted that a search for matching 
circles has been performed in Ref. [ll| and Ref. 1^. 
The authors found that their circle matching statistic 
was degraded considerably by the ISW and velocity con- 
tributions to the temperature fluctuations in the CMB. 
In Ref. [ij, it was found that the matching circles of 
a Poincare dodecahedron universe cannot be found by 
current CMB searches because of this degradation. In 
Ref. [i3, the search included filtering and combines sets 
of expected circles, but disregarded the m weighting in 
the S statistic, and provides no false positive thresholds. 
As we have shown, the use of power spectrum filtering, 
the combination search, and the proper S statistic that 
includes m weighting can indeed probe beyond contami- 
nating factors to determine if there are matching circles 
in the CMB. Applying these techniques to the putative 
Poincare Dodecahedral model produced no statistically 
significant sets of matching circles. Since the smallest 
matching circles in our simulated sky have radii of 16.4°, 
we have not been able to directly establish a lower limit 
on the angular radii that can be probed by our combi- 
nation search. However, we can extrapolate from the 
results shown in Figures and |S1 to put the lower limit 
at a ~ 5° . This effectively rules out the Poincare Dodec- 
ahedral model as an interesting shape for the Universe. 



XI. CONCLUSION 

The claim that the topology of our Universe had been 
found to be that of the Poincare dodecahedron does not 
stand up under scrutiny. The signature found in Ref. |gj 
disappears when one uses the proper S statistic and con- 



siders the false positive threshold. While the shape of our 
Universe remains a mystery, the matching circles test can 
be used to place a lower bound on the size of the Universe. 
The previous limit of 24 Gpc ^ can be extended by about 
10% using filtering of the WMAP power spectrum. A full 
search with optimal filtering is now underway. 
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